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Abstract 

We describe a basic framework for studying dynamic scaling that 
has roots in dynamical systems and probability theory. Within this 
framework, we study Smoluchowski's coagulation equation for the three 
simplest rate kernels K{x,y) — 2, x + y and xy. In another work, we 
classified all self-similar solutions and all universality classes (domains 
of attraction) for scaling limits under weak convergence (Comm. Pure 
Appl. Math 57 (2004)1197-1232). Here we add to this a complete de- 
scription of the set of all limit points of solutions modulo scaling (the 
scaling attractor) and the dynamics on this limit set (the ultimate dy- 
namics). The main tool is Bertoin's Levy-Khintchine representation 
formula for eternal solutions of Smoluchowski's equation (Adv. Appl. 
Prob. 12 (2002) 547-64). This representation linearizes the dynamics 
on the scaling attractor, revealing these dynamics to be conjugate to a 
continuous dilation, and chaotic in a classical sense. Furthermore, our 
study of scaling limits explains how Smoluchowski dynamics "compact- 
ifies" in a natural way that accounts for clusters of zero and infinite 
size (dust and gel). 
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1 Introduction 

Smoluchowski's coagulation equation is a fundamental mean-field model of 
clustering processes. The merging of clusters of mass x and mass y to 
produce clusters of mass x + y occurs at a mass-action rate modulated by 
a symmetric rate kernel K{x,y). Formally, the evolution equation for the 
density n{t, x) of the size distribution reads 



Many kernels arising in applications are homogeneous, that is, there is 7 such 
that K(ax,ay) = a'^K{x,y) for every a,x,y > 0. Wc restrict attention to 
the "solvable" kernels K{x, y) = 2, x + y and xy (7 = 0,1 and 2 respectively) 
which may be studied via the Laplace transform. These special kernels arise 
in a variety of applications, including the aggregation of colloids [6] (A' = 2), 
droplet formation in clouds [11], a phase transition for parking [7], Burgers' 
model of turbulence [3] {K = x + y), random graphs [1] {K = xy), and 
kinetics of polymerization [18] (all kernels). 

As time proceeds, the typical cluster size grows, and an issue of relevance 
for homogeneous kernels is whether and how the size distribution develops 
toward self-similar form. In the physics literature, this is called the dynamic 
scaling problem. This work continues our study of dynamic scaling for 
Smoluchowski's equation for the solvable kernels. We lay out a general 
framework for the analysis of dynamic scaling that is inspired by elements 
from both dynamical systems and probability theory. The main issues may 
be set out list of basic and general questions: 

• What sealing solutions exist? Here we seek self-similar solutions, or 
fixed points of the dynamics modulo scaling. 

• What are the domains of attraction of these scaling solutions? These 

comprise the universality classes for dynamic sealing. 

• What limit points are possible under scaling dynamics in general? We 
call the set of such points the scaling attractor of the system. 

• How can we describe the dynamics on the scaling attractor? We call 
this the ultimate dynamics of the system. 



dtn{t,x) = - K{x -y,y)n{t,x -y)n{t,y)dy 




(1.1) 



• How complicated can the ultimate dynamics be? 



G. Menon and R. L. Pego 



3 



While evidently stated in dynamical terms, strong motivation for this 
framework comes from the classical limit theorems of probability theory 
developed by the pioneers of the subject in the 1930s. These theorems 
concern limits of scaled sums of independent and identically distributed 
random variables Xi, X2, . . . with common distribution F. The distribution 
of = X]"=i is the n-fold convolution of F with itself and this can be 
regarded as a discrete evolution law for F. The theory as laid out in the 
marvelous book of Feller [9] provides complete answers to the questions 
above: 

• Scaling solutions. The normal distribution is the unique invariant dis- 
tribution of finite variance. More generally, all invariant distributions 
(including those with heavy tails) comprise a two-parameter family, 
the stable laws first characterized completely by Levy (see [14] for a 
historical account). 

• Domains of attraction. The central limit theorem states that the nor- 
mal law attracts all distributions with finite variance. The domains of 
attraction of the stable laws are completely classified in terms of the 
power-law behavior (more precisely, regular variation) and skewness 
of the 2nd-moment distribution function. There are no other invariant 
distributions or domains of attraction in the limit n — > 00. 

• Infinite divisibility. The most general distributions that can arise as 
rescaled limits for some subsequence Uj — 00 are the infinitely divisible 
distributions, characterized by the famous Levy-Khintchine formula. 
The characteristic function (p{k) = K{e'^^^) of an infinitely divisible 
random variable X is of the form <p{k) = exp^'(fc) where 



with jjj(l A X ^)M(cix) < 00 and 6 G M. The measure M is called the 
canonical measure. 

• Ultimate dynamics. Under appropriate rescaling, nondegenerate limit 
points for the discrete evolution n 1-^ Sn correspond one-to-one to 
Levy processes (continuous-time random walks with stationary inde- 
pendent increments). These are stochastic processes Xf {t > 0) with 
right continuous paths with left limits, obeying the semigroup formula 
(fitik) = E{e^'''^*) = e**^*^). Rescaling amplitude and t (scaling dynam- 
ics) corresponds to linear tranformations (dilation, stretching, shifts) 
of the canonical measure. 




,ikx 
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• Chaos. There exist distributions {Doeblin's universal laws) for which 
every infinitely divisible distribution appears as some limit point for 
rescaled self-convolutions. 

In an earlier article [16], wc characterized the approach to self-similarity 
for Smoluchowski's equation as t ^ oo for K = 2 and x+y, and the approach 
to self-similar blow-up as t approaches the gelation time for K = xy. By 
way of addressing the first two issues in the framework above we found: 

• For each solvable kernel, with degree of homogeneity 7 G {0,1,2}, 
there is a (classically known) self-similar solution with finite 7 -|- 1st 
moment, unique up to normalization. This solution has exponential 
decay as x ^ 00. But more generally, there is a one-parameter family 
of self-similar solutions, corresponding to distribution functions writ- 
ten Fp^^ for p G (0, 1], 7 G {0, 1, 2}. We recover the classical self-similar 
solution at the endpoint p = 1. For < p < 1, the Fp,7 have infinite 
7 -|- 1st moment and are directly related to important heavy-tailed 
distributions of probability theory — Mittag-Lefiler distributions for 
K = 2, and stable laws of maximum skewness ior K = x + y and xy. 
For K = X + y these solutions were first discovered by Bertoin by a 
different argument [4]. 

• The classical self-similar solution with p = 1 attracts all solutions 
with finite 7 + 1st moment. In general, the domains of attraction of 
self-similar solutions are characterized by the regular variation of the 
7-I- 1st moment distribution. A positive measure vq lies in the domain 
of attraction of the self-similar solution Fp^^ if and only if 



for some function L slowly varying at 00 {L{\x) / L{\) — > 1 as A — > 00 
for all X > 0). There are no other self-similar solutions or domains of 
attraction. 

Our present goal is to investigate more generally the scaling dynamics 
of Smoluchowski's equation for the solvable kernels, as t ^ 00 for K = 2 
and x + y and as t approaches the gelation time for K = xy. To stress 
the probabilistic analogy, if our earlier article was a study of stable laws, 
here we study infinite divisibility. The basis for our work is Bertoin's char- 
acterization of eternal solutions for Smoluchowski's equation with additive 
kernel K = x + y [4]. Eternal solutions for this kernel are defined for all 




X 



00 



(1.2) 
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t G (—00, co) (i.e., they may be extended backwards in time globally), and 
thus they are 'infinitely divisible' under clustering dynamics. Bertoin estab- 
lished a remarkable Levy-Khintchine-type representation for these solutions. 
We generalize this result to other solvable kernels, and revisit it from a dy- 
namical systems viewpoint. In the context of the framework above, we find 
that 

• The (proper) scaling attractor corresponds in one-to-one fashion with 
eternal solutions of Smoluchowski's equation, and these have a Levy- 
Khintchine representation for each solvable kernel. 

• The Levy-Khintchine representation linearizes the dynamics on the 
attractor. Nonlinear evolution by Smoluchowski's equation on the at- 
tractor is conjugate to a group of simple scaling transformations on 
the measures that generate the representation. Heuristically, the clas- 
sification of the domains of attraction in [16] shows that the dynamics 
are extremely sensitive to the tails of the initial size distribution. The 
Levy-Khintchine representation makes this precise: the ultimate dy- 
namics on the scaling attractor is conjugate to a continuous dilation 
map. 

• One may use the Levy-Khintchine representation to construct orbits 
with complicated dynamics. The scaling attractor contains a dense 
family of scaling-periodic solutions. Furthermore, there are eternal 
solutions with trajectories dense in the scaling attractor — we call 
these Doeblin solutions. And, for any given scaling trajectory, there is 
a dense set of initial data whose forward trajectories shadow the given 
one. 

In addition, this study of scaling limits reveals how Smoluchowski dy- 
namics "compactifics" in a natural way that accounts for clusters of zero and 
infinite size (dust and gel). Considering defective limits on (0, 00) that con- 
centrate probability at and 00 yields a well-posed dynamics of "extended 
solutions" on [O. oo]. Proper solutions remain fundamental, but considering 
extended solutions with dust and gel helps to understand just how the tails 
of initial data determine long-time behavior. 

We remark that scaling-periodic solutions are analogous to the semi- 
stable laws in probability theory [15]. Our "Doeblin solutions" are con- 
structed by "packing the tails" of the Levy measure in a fashion entirely 
analogous to the construction of Doeblin's universal laws in probability [9, 
XVII. 9]. In this connection, it is interesting to note that the examples 
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in [9, XVII.9], dismissed by Feller as "primarily of curiosity value," closely 
resemble modern treatments of chaos, and Doeblin's construction appears 
particularly prescient. 

The solvable cases of Smoluchowski's equation correspond to sophisti- 
cated stochastic models with a rich theory (see [1, 5] for excellent reviews), 
so perhaps it is no accident that the classical probabilistic methods work 
so well. But let us stress that our work really relies only on the analytical 
methods for studying scaling limits that lie behind the classical limit theo- 
rems. These methods are simple and powerful and should be of utility for 
understanding scaling phenomena in other applications that have no obvious 
probabilistic meaning. Thus, no knowledge of probability is presumed and 
(almost) all details are included (though there is no substitute for reading 
Feller!). 



2 Statement of results 

In this section, we state our results precisely in a setting that unifies the 
treatment of dynamic scaling for all the solvable kernels. 

Let E denote the open interval (0, oo), M the space of non- negative 
Radon measures on E, and V the space of probability measures on E. We 
will always use the weak topology on M. and V. We also let E denote the 
closed half line with point at infinity, E = [0, oo] = [0, oo) U {oo}, and let P 
be the space of probability measures on E. 

Rigorous theories for solutions where ut{dx) = n{t,x)dx is a general 
size-distribution measure on E = (0, oo), thus accounting for both continu- 
ous and discrete size distributions in one general setting, are based on the 
moment identity 

^Jj{x)ut{dx)= ^j^jj{x,y)K{x,y)ut{dx)vt{dy), (2.1) 

where / is a suitable test function and /(x, y) = f{x + y) — f{x) — f{y), see 
[16, 17, 10]. Let mg{t) := j^x^ut{dx) denote the 0-ik moment of ut. By the 
results of [16], for a solvable kernel of homogeneity 7, any initial measure 
uto with finite 7-th moment m-y(io) determines a unique continuous weak 
solution 

v={vt{dx),te[tG,T^^)). (2.2) 
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For convenience we can always scale the initial data so that 
' I {K = 2), 

to := < {K = x + y), and m-y(to) = / x^vtoidx) = 1. (2.3) 
^ -1 {K = xy), 

Then Tmax = cxd for = 2 and x + y, and Tmax = for K = xy. For each 
solvable kernel, m-y{t) = y'^vt{dy) is an explicitly known function — from 
(2.1) with f{x) = we find 



r r 



{K = 2), 



m^{t) = < 1 (iC = X + y), 
^ \t\-^ {K = xy). 



(2.4) 



The solution ut is typically not a probability measure because the total 
number of clusters decreases in time, but there is a naturally associated 
probability measure Ft{dx) with distribution function Ft(x) defined by 



Ft{x) = [ y-^vtidy) I [ y-'i^tidy). 

J{0,x] I JE 



(2.5) 



In this way, we regard Smoluchowski's equation as defining a continuous 
dynamical system on the phase space V. 



2.1 Eternal solutions 

For exceptional initial data we may also solve backwards in time (mean- 
ing is divisible under clustering dynamics). The maximum possible in- 
terval of existence that can be obtained in this way is denoted (Tmin, Trnax); 
where Tmin,7max depend only on the kernel and x'^Vfoidx). With the 
normalization (2.3), the maximum possible interval of existence turns out 
to be 

r (0,oo) {K = 2), 

(rmin, Tmax) = < (-00, Oo) {K = X + y), (2.6) 

[ (-00,0) {K = xy). 

Solutions which are defined on this maximum interval of existence are the 
analog of infinitely divisible laws in probability. 

Definition 2.1. A solution to Smoluchowski's equation that is defined for 
all t G (Tmin, Tmax) IS Called an eternal solution. 
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2.2 The scaling attractor 

A central idea in dynamical systems theory is to understand the long-time 
behavior of solutions through the notions of attractor and w-limit sets. Coag- 
ulation equations transport mass from small to large scales, and all mass es- 
capes as i — > Tmax- To obtain non-trivial long-time behavior wc must rescale 
solutions. We adopt the following definitions for such scaling dynamical sys- 
tems. Below, Tn G [toi^max), Pn > 0. We will often use the same letter to 
denote a measure and its distribution function, e.g., F{x) = Jj^^j F{dx). 

Definition 2.2. The (proper) scaling u-limit set of a solution u to Smolu- 

chowski's equation is the set of probability measures F on E for which there 
exist sequences T„ Tmax, Pn oo, such that FT„{Pnx) — > F{x) at every 
point of continuity of F. 

Definition 2.3. The (proper) scaling attractor, Ap, is the set of proba- 
bility measures F on E for which there exists a sequence of solutions i^^"^ 
defined for t G [toj^max), and sequences — > T'maxj Pn — ^ OO; such that 
Fj^{Pnx) — > F{x) at every point of continuity of F. 

As a consequence of continuous dependence of solutions on initial data 
(forward and backward in time), we will show that the scaling attractor is 
an invariant set, and that points on the proper scaling attractor and eternal 
solutions are in one-to-one correspondence. 

Theorem 2.1. (a) The proper scaling attractor Ap is invariant: If v is 
a solution of Smoluchowski 's equation, and Ft G Ap for some t, then 
the solution is eternal and Ft G Ap for all t G (Tmin, ^max)- 

(b) A probability measure F belongs to Ap if and only if F{dx) = x'^vt^{dx) 
for some eternal solution v, where to is as in (2.3). 

The perfect definition of an attractor remains elusive (see for example, 
the discussion in [12, Ch. 1.6]). Definition 2.3 is perhaps the simplest for 
dynamical systems. It also has the virtue of generalizing the probabilistic 
notion of domains of partial attraction [9, XVII]. However, some typical 
properties that hold in finite-dimensional dynamical systems do not hold 
here. For example, it need not be the case that every solution has a non- 
empty scaling w-limit set. Nor is closed. Defective limits are possible, as 
shown in [16]. See [13] for a discussion of related issues in the probabilistic 
context. It turns out that we can cure these defects and account for limits 
that involve mass concentrating at zero or leaking to infinity, by the simple 
expedient of allowing limits to be probability measures on £^ = [0, oo] . 



G. Menon and R. L. Pego 



9 



Definition 2.4. The full scaling oj-limit set of a solution v to Smolu- 

chowski's equation is the set of probability measures F on E with the prop- 
erty in Definition 2.2. The full scaling attractor, A, is the set of probability 
measures F on E with the property in Definition 2.3. 

The space V of probability measures on E, equipped with the weak 
topology, is compact — any sequence contains a converging subsequence. 
We will show that Smoluchowski dynamics naturally extends by continuity 
from V to V. Such "extended solutions" have probability distributions that 
may include atoms at and oo, allowing for the possibility that clusters 
have zero size ("dust") or infinite size ("gel") with positive probability. To 
interpret these physically, one should recognize that of course and oo are 
idealizations relative to a given scale of measuring cluster size. 

We defer detailed discussion of extended solutions to section 9. There 
we extend Theorem 2.1 to relate the full scaling attractor A (now a compact 
set that is the closure of .4p) to the set of eternal extended solutions. Also, 
the Levy-Khintchine representation and linearization theorems in the next 
two sections have elegant extensions involving extended solutions. First, 
however, we think it appropriate to focus on standard weak solutions, and 
develop the theory without dust in our eyes, so to speak. 

2.3 Levy-Khintchine representations 

In probability theory, infinitely divisible distributions are parametrized by 
the Levy-Khintchine representation theorem, which expresses the log of the 
characteristic function (Fourier transform) in terms of a measure that sat- 
isfies certain finiteness conditions. In particular [9, XIII. 7], a function uj{q) 
is the Laplace transform Jq e~'^^ F{dx) of an infinitely divisible probability 
measure F supported on [0, oo) if and only if io{q) = exp(— $(^)) where the 
Laplace exponent $ admits the representation 




(2.7) 



for some measure G on [0, oo) that satisfies 




and 




for all a; > 0. 



(2.8) 



Equivalently, 




(2.9) 



G. Menon and R. L. Pego 



10 



We need a name for measures with this property, although none seems stan- 
dard. Such measures determine the jump-size distribTition for subordinators 
— increasing continuous-time random walks with stationary independent 
increments — hence we call them s-measures. To handle defective limits, it 
is convenient to allow y~^G{dy) to have an atom at oo. 

Definition 2.5. A measure G on [0, oo) is an s-measure if (2.9) holds. A 
pair (G, 5oo) is called an s-measure on E = [0, oo] if G is an s-measure and 

5oo > 0. goo is called the charge of y~^G{dy) at oo, and we will abuse 
notation by denoting the pair {G,goo) by G. In addition, we say that an 
s-measure (or s-measure) G is divergent if 

G(0) > or / y^^G{dy) = oo. (2.10) 
Je 

Here recall we use the notation G{x) = G{dy). If g^o = wc identify 
G with (G, 0). The space of s-measures has a natural weak topology which 
will prove fundamental in our study of scaling dynamics. 

Definition 2.6. A sequence of s-measures converges to an s-measure 
G as n — > oo, if at every point x G (0, oo) of continuity of G we have 



/ 

J\o. 



G(-){dy) [ G{dy) (2.11) 



and 

/ y-^G(")(dy) ^ / y-^G{dy). (2.12) 

J[x,oo\ J{x,(xi\ 

The integrals in (2.12) include the charge at oo, if any. We note that 
in view of the weak convergence implied by (2.11), convergence of to 
an s-measure G (having g^ = 0) is equivalent to (2.11) together with the 
tightness condition 

/ y~^G^'^\dy) as x ^ oo, uniformly in n. (2-13) 

J [XjOo] 

Bertoin's main theorem in [4] shows that eternal solutions for K = x + y 
are in one-to-one correspondence with divergent s-measures. (More precisely, 
Bertoin formulates his result in terms of "Levy pairs," separating the atom 
at the origin from a jump measure on (0, oo)). We extend this result as 
follows. Let ly be an arbitrary solution to Smoluchowski's equation for a 
solvable kernel of homogeneity 7. Since the 7-th moment of Vf is finite. 
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x'^'^^i>t{dx) is an s-measure. Rescaling, we associate with ut the s-measure 
Gt defined by 



Gt{dx) = x^+^ut{X{t)dx), X{t) 



1 iK = 2), 

e* {K = x + y), (2.14) 

{K = xy). 



Our choice of rescaUng ensures that if the total measure Gt{E) is finite for 
some t, then it is constant: Gt{E) = m^+i(i)/A(i)'*'+"^ =const. 
One computes that 

/ y-^Gt{dy) = '^ = { C-* {K = x + y), (2.15) 
•^0 -^W^ [ \t\ {K = xy) 

The well-posedness theorem [16] impHes that solutions of Smoluchowski's 
equation normalized according to (2.3) that exist on any time interval [t, ^?max) 
are in one-to-one correspondence with s-measures G that satisfy 

/•oo 

G(0) = and / y-^G{dy) = my{t)/X{ty , 
Jo 

via G = Gt- Through studying the limit t [ Tmin, we find that eternal 
solutions may be characterized as follows. 

Theorem 2.2. (a) Let v he an eternal solution for Smoluchowski's equa- 
tion with K = 2,x -\- y or xy. Then there is a divergent s-measure H 
such that Gt converges to H as t [T^ 



'-mm- 



(b) Conversely, for every divergent s-measure H , there is a unique eternal 
solution V such that Gt converges to H as t [ Tj^in- 

(c) Let ©pi Ap Sd map the (proper) scaling attractor Ap to the set 
Sd of divergent s-measures by (5p{F) = H, where H is the diver- 
gent s-measure associated to the eternal solution v such that F{dx) = 
x'^L'toidx) with to as in (2.3). Then &p is a bi-continuous bijection. 

The procedure for obtaining the eternal solution i/ from the divergent s- 
measure H is nonlinear and is different for each kernel (see Theorems 4.3, 5.4, 
and 6.1 below). It seems natural to call Theorem 2.2 a Levy-Khintchine 
representation for the scaling attractor Ap — as we will see, eternal solutions 
are determined through the Laplace exponents associated with divergent s- 
measures. 

In section 9, the correspondence in Theorem 2.2 is expanded to one 
between eternal extended solutions and arbitrary s-measures. 
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2.4 Linearization of ultimate dynamics 

There are two natural group actions on the class of eternal solutions that are 
related to scaling dynamics, arising from time evolution and rescaling of size. 
A straightforward but remarkable consequence of the scaling properties of 
Smoluchowski's equation is that nonlinear dynamics (time evolution) on the 
scaling attractor Ap is conjugate to a simple linear scaling transformation 
of the divergent s-measures that correspond by Theorem 2.2. 

Theorem 2.3. Let u be a solution of Smoluchowski's equation with K = 2, 
X + y or xy. Given scaling parameters a > and b > 0, let 

' auatibdx) iK = 2), 
utidx) = < bvt+iogaibdx) {K = x + y), (2.16) 
^ ab'^Uat{bdx) {K = xy), 

with associated probability distribution function 

j Patibx) {K = 2 or xy), . . 

^^^'^^-X Fi+iog,(6x) {K = x + y). ^'■'^> 

Then v is again a solution. If v is eternal and H its associated divergent 
s-measure, then u is eternal and its associated divergent s-measure is given 
by 

( ab-^H{bx) {K = 2), 

H{x) = I aH-^H{a-^bx) [K = x + y), (2.18) 
[ a-H-^H{abx) {K = xy). 

Proof. The proof is simple, based on Theorem 2.2 and the scaling properties 
of Smoluchowski's equation. First, one checks that (2.16) determines a solu- 
tion, by scaling the moment identity (2.1) in each case. Next, compute that 
if the s-measure Gt is associated with vt as in (2.14), then the corresponding 
s-measure associated with i>j is given by 

{xvt{dx) = ab-^Gatibdx) iK = 2), 

x'^i>t{e* dx) = a^b~^Gt+ioga{a~^bdx) {K = x + y), 
x^i>t{\t\-^dx) = a-H-^Gat{abdx) {K = xy). 

(2.19) 

Then take t j Tmm and apply Theorem 2.2 to deduce (2.18). □ 



Theorem 2.4. Let u be an eternal solution with corresponding divergent 
s-measure H and let Ft be as in (2.5) for K = 2, x + y or xy. For each 



G. Menon and R. L. Pego 



13 



t £ {TmimT^ax), let Hf = &p{Ft) be the divergent s-measure associated to 
Ft^A^. Then 

( tH{x) {K = 2), 

Ht{x) = I e^^H{e-^x) {K = x + y), (2.20) 
( \t\-^H{\t\x) {K = xy). 

Proof. Take 6=1 and put t = to in (2.17), then substitute a = t, e*, \t\ for 
K = 2, x + y, xy respectively to obtain Ft^ = Ft- Then the corresponding 
divergent s-measure H = &p{Ftf^) is found from (2.18). □ 

By this theorem, we see that in terms of the divergent s-measure that 
corresponds to the solution, the time evolution on the scaling attractor Ap 
is governed by the linear equations 

tdtHt = Ht {K = 2), (2.21) 
idt + xd,)Ht = 2Ht {K = x + y), (2.22) 
{tdt - xd^)Ht = -2Ht {K = xy). (2.23) 

2.5 How initial tails encode scaling limits 

The long-time scaling behavior is very sensitive to the initial distribution 
of the largest clusters in the system, as indicated by the characterization 
of domains of attraction via (1.2), and the linearization theorem above. In 
fact, the long-time scaling dynamics is encoded in the tails of initial data in 
a simple fashion related to the Levy-Khintchine representation. 

Theorem 2.5. Let F £ Ap with associated divergent s-measure H . Let 
i/^"^ he any sequence of solutions defined for t > to, with associated initial 
s-measures given by G^'^\dx) = x'^'^^i'^^\dx). Let — > Tmax; l^n oo. 
Then the following are equivalent: 

(i) Fj^{(3nx) — > F{x) as n ^ oo, at every point of continuity. 

(ii) The rescaled initial s-measures defined by 

' p-'TnG(^\Pnx) {K = 2), 

G(")(x) = i /3-ie2^"G('^)(e-^"/3„x) {K = x-{-y), (2.24) 
^ p-'\Tn\-^G(-\\Tn\Pnx) {K = xy), 

have the property that G^"^ converges to H as n ^ oo. 
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This result generalizes to the fuh attractor A, with H replaced by the corre- 
sponding s-measure, see section 9. We remark that in the proof it is shown 
that for the convergence in part (ii) to hold, it is necessary that e~^"/?„ — > oo 
ioi K = X + y, and |T„|/3„ — >■ oo for = xy. 

2.6 Signatures of chaos 

The dilational representation of dynamics in (2.20) in terms of the Levy- 
Khintchine representation means that Smoluchowski dynamics on the scal- 
ing attractor is a continuous analog of a Bernoulh shift map, a classical 
paradigm for chaotic dynamics. We demonstrate the utility of this repre- 
sentation by constructing solutions with both chaotic and regular orbits, 
and by proving a shadowing theorem illustrating sensitive dependence on 
the tails. 

2.6.1 Solutions with dense limit sets 

Theorem 2.6. There exists an eternal solution v whose scaling ui-limit set 
contains every element of the full scaling attractor A- 

We call such solutions Doeblin solutions by analogy with Doeblin's univer- 
sal laws. The construction follows Feller closely and relies only on general 
principles (separability of V and S^, and continuity of the bijection Sp). 
Theorem 2.6 tells us that A cannot be decomposed into invariant subsets. 

2.6.2 Scaling periodic solutions 

Another classical signature of chaos is the density of periodic solutions. The 
notion of periodicity generalizes as follows. 

Definition 2.7. Let v he a solution and define Ft{x) by (2.5). We say ly is 
scaling-periodic if for some ti > to and P > 1, 

Ft, {/3x) = Ft,, {x) for all x > 0. (2.25) 

These are analogous to semi-stable laws in probability theory [15]. The 
Levy-Khintchine representation yields a simple characterization. 

Theorem 2.7. A scaling-periodic solution of Smoluchowski' s equation with 
kernel K = 2,x-\-y or xy is eternal, and its divergent s-measure H satisfies 

H{x) = aH{bx) (2.26) 

for some a > 0, 6 > 1 such that either 
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(i) a = 1 and H is an atom at the origin, or 

(a) a < 1 and ab > 1. 

Conversely, if H is a measure on [0, oo) with H{x) = aH{hx), where a > 
0, 6 > 1 and (i) or (ii) hold, then H is a divergent s-measure and the 
corresponding eternal solution is scaling-periodic. 

Case (i) is simple but important. The corresponding scaling-periodic 
solutions are the self-similar solutions with exponential decay. More gen- 
erally, all self-similar solutions are determined by divergent s-measures of 
the power-law form H{x) = CpX^~P, < p < 1. Scaling-periodic solutions 
that are not self-similar solutions are generated by (ii). Thus, there are un- 
countably many scaling-periodic solutions. Moreover, the Levy-Khintchine 
representation allows us to prove 

Theorem 2.8. Scaling-periodic solutions are dense in A. 



2.6.3 Shadowing and sensitive dependence on the initial tails 

We show that asymptotically similar initial tails imply shadowing of scaled 
solution trajectories. To study shadowing, we note that the space V of 
probability measures on E is metrizable and compact. We let dist(-, •) denote 
any metric on V which induces the weak topology. 

Theorem 2.9. Let v and v denote any two solutions of Smoluchowski's 
equation defined on [ioj^max); o,nd let the associated initial s-measures be 

G{dx) = x'^+^ut^idx), G{dx) = x^+^vt^x), (2.27) 

with Laplace exponents and (p respectively associated as in (2.7). Assume 
that 

^(g)Mg) ~ L(l/g) asq^O, (2.28) 

where L is slowly varying at oo. Suppose that b{t) t oo as t j Tj^ax, o,nd 
define 

' {t/L{b), b) {K = 2), 

{t,b) = l {t- log L(6e-*), b/L{be-*) ) {K = x + y), (2.29) 
. {tL{\t\b), bLi\t\b) ) {K = xy), 

so that b/X(t) = b/X{t) with X{t) as in (2.14). Then we have 

dist{Ft{b{t)dx),Ft{b{t)dx)) ^0 as t ^ T^e^. (2.30) 
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The simplest situation, requiring no readjustment of the scahng {i = t, 
b = b) is when L = 1. When condition (2.28) holds, the solutions u and i' 
have identical scaling w-limit sets, for example. If one of the solutions in 
this theorem, say v', is self-similar, then the sufficient condition (2.28) for 
shadowing in this theorem is equivalent to (1.2) (see [16], (5.3) and (5.7) for 
K = 2, and (7.2) and (7.4) fox K = x + y). Hence (2.28) is also necessary, 
according to the classification theorem on domains of attraction. It appears 
that in general the sufficient condition for shadowing given in this theorem 
may not always be necessary. But we will not pursue this issue here. 

The sensitivity of solutions to initial tails in the weak topology is revealed 
strikingly in Theorem 2.9. The topology of weak convergence is undoubtedly 
natural for limit theorems, for example the approach to self-similarity. On 
the other hand, this topology cannot distinguish the tails, as the following 
"cut-and-pastc" argument shows. Let F = x'^uto and F = x'^vt^ be given 
initial data for two solutions, and define 

Then ^ F as ra ^ 00, and F^'^^ has Laplace exponent given by 

. /" (1^) + r (i^) ,p^,y, (2.32) 

q Jo \ Qy J Jn \ Qy J 

from which one sees easily that if G{E) = oo, then Lp^^\q) ~ (^((7) as g ^ 0. 
Thus, according to Theorem 2.9, the solution z^^"^ generated by F'^"-' will 
shadow u. This justifies the statement made in the introduction that for any 
given scaling trajectory, there is a dense set of initial data whose forward 
trajectories shadow the given one. 

2.7 Plan of the paper 

In section 3 we establish some basic facts regarding convergence of the mea- 
sures that generate the Levy-Khintchine representation. The analysis of 
eternal solutions is different in detail for the constant and additive kernels, 
so we treat these cases in turn, establishing Theorems 2.1, 2.2, and 2.5 for 
these kernels in sections 4 and 5. The multiplicative case reduces mathe- 
matically to the additive by a change of variables, and is treated briefly in 
section 6. We emphasize that the results in this case concern the behavior of 
solutions approaching the gelation time, so perhaps this is the case of most 
interest physically. 
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With the Levy-Khintchine representation in hand, we then construct 

Doebhn solutions in section 7 and scaling-pcriodic solutions in section 8. 
Extended solutions and the full scaling attractor A arc studied in section 
9, and the shadowing theorem 2.9 is proved in section 10, where we also 
provide a streamlined treatment of the domains of attraction. 

3 Laplace exponents and limits of s-measures 

The main analytic tool in the study of the solvable kernels is the Laplace 
transform. Recall that a sequence of probability measures F^") is said to 
converge weakly to a probability measure F if the distribution functions 
— > F{x) at every point of continuity of the limit. It is basic [9, 
XIII. 1] thatFW converges weakly to F if and only if the Laplace transforms 
converge pointwise: 

/ e-9^F('^)(dx) ^ / e-«^F(dx), for all g > 0. 
Je Je 

We will need the following refinements of this result for s-measures. With 
any s-measure G we associate "Laplace exponents" $ and * (with $ = dq^) 
defined for g G C+ = {A G C : ReA > 0} by 

m = / G{dx), ^{q) = / ^±^G{dx). (3.1) 

Je ^ JE ^ 

If go and Qoo denote amplitudes of the atoms of the measure (1 A y~^)G{dy) 
at and oo, respectively, this means 

r 1 _ Q-i^ 

^{q) = q9o + 9oc+ G{dx), (3.2) 

J(0,oo) X 

*(«) = 7^q'9o + qgoo+ -^G{dx). (3.3) 

We use the terminology Laplace exponent in accordance with probabilists' 
usage. If we need to distinguish the two types of exponents, we will refer 
to <I> and ^' as Laplace exponents of the first and second order, respectively. 
Observe that 

dg^ = d'^^= [ e-«^G((ix). (3.4) 

J[0,oo) 

These functions are Laplace transforms of a positive measure, thus are com- 
pletely monotone functions on (0, oo). 
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We note that the amphtude of the atom of (1 A y ^)G{dy) at oo is 

Qoo = Hm ^{q) = hm (3.5) 
g— >0+ q— >0+ 

thus the s-measure G is an s-measure if and only if this vanishes. Further- 
more, we claim that G is divergent if and only if 

lim $(g) = cxD, equivalently lim q~^'^{q) = oo. (3.6) 

g— >oo q—>co 

To prove this, observe that 

Hq) <q9o+ f x-^G{dx). 

J(0,oo) 

Thus, if lim^^oo ^{q) = oo, then G satisfies (2.10). Conversely, if G satisfies 
(2.10), then limg^oo ^{q) = c)o by the monotone convergence theorem. The 
proof for ^ is similar. We integrate by parts and use Fubini's theorem to 
obtain 



q-'^{q) = iqgo+g^+ [ (1 - g-^^) / y-^G{dy)dx 
^ J{0,oo) J{x,oo) 

< \q9Q + 9oo+ I y~^G{dy). 

^ J(0,oo) 



'(0,oo) 

Thus, if q~^^{q) oo then G satisfies (2.10). The converse follows from 
the monotone convergence theorem. 

Theorem 3.1. Let G'^"-' be a sequence of's-measures with Laplace exponents 
$H and*H. Then, taking n — > oo, the following are equivalent: 

fi) converges to an s-measure G with Laplace exponents $ and 

(ii) ^{q) := lim^^oo ^^"\q) exists for each q > 0. 

(Hi) "^(q) ■= lim„_>oo ^^'^\q) exists for each q > 0. 

Proof, (i) implies (ii): Fix q> and let e > 0. (2.12) allows us to choose a 
such that o is a point of continuity for G and for every n 

[ e-«^x-^ (G(") (dx) + G{dx)) < e-««C < e. 

J [a,oo] 

On the other hand, (2.11) guarantees 

/V - e"^'')x'^G'^''\dx) ^ /V - e"^'')x-^G{dx). 
Jo Jo 
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Using again (2.12), we conclude that for large n, \^^'"\q) — ^{q)\ < £. 

(ii) implies (i): 1. Claim: <5 is analytic in C+ and — > $ uniformly 
on compact subsets of C-|-. Proof: Let K C C+ be compact. The claim 
follows from the estimate: 

supsup|$(")(g)| < oo. (3.7) 

n qeK 

Indeed, by Montcl's theorem, (3.7) implies {^^'"'^}'^^-^ are a normal family of 
analytic functions (i.e., precompact in the uniform topology). Thus, every 
subsequence has a further subsequence converging uniformly to an analytic 
function. Since every subsequence converges (pointwise) to this implies 
(|)(") $ uniformly and <I> is analytic. It remains to prove (3.7). We 
integrate by parts to obtain 

g-i$(")(5) = g(")(0) + / e""^ ( / 2/-^G(")(dy) I dx. 

Je J 

Thus, for any a > 0, supReg>„ $(")(?) | < a-i$(")(a). Since $(")(g) 
converges for all q > 0, we have sup„ supRgq>„ |g~^$'^")(g)| < oo. This 
proves (3.7). 

2. Cauchy's integral formula and the claim imply d^^^"'^ — >■ for 
every € N. Since are completely monotone, so is the limit dq^. 

3. Thus, dq^ = JjQ e~'^^G{dx) is the Laplace transform of a measure 

G on [0, oo). We integrate with respect to q with (/oo := '^(0^) defined to 
be the charge of y~^G{dy) at oo, and use Tonclli's theorem to obtain (3.1). 
Note that Jj^ x~^G{dx) < oo because ^{q) < oo for each fixed q, so G is 
an s-measure. 

4- The convergence 9^$^") — > dq^ is equivalent to weak convergence of 
G*-"-* to G on [0, oo), meaning (2.11) holds. This implies that for every point 
X of continuity of G, as n ^ oo we have 

/ (l-e-'^y)y-'G^^\dy)-^ [ {1 - e-'^y)y-'G{dy), (3.8) 

J[0,x] J[0,x] 

and together with (ii) this yields that Jj^ y~^G^'^\dy) is bounded and 

[ e-iyy-^G^''\dy) [ e-'iyy-^G{dy). 

J[x,oo\ J[x,oo\ 

From (ii) then follows (2.12). This proves G^") G. 
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(ii) implies (iii): This is due to = $^"^(5) ds and monotonicity. 

(iii) implies (ii): Since *(")(g) = J^^^''\s)ds > y^^''\lq), we find 
that (3.7) holds as in step 1 above. Then for every subsequence of ^^^^ 
there is a further subsequence that converges on compact sets of to an 
analytic limit This limit is unique due to (iii), and (ii) follows. (It follows 
also that ^(") ^ uniformly on compact sets.) □ 



4 The constant kernel 

In this section we study eternal solutions and the Levy-Khintchine represen- 
tation in particular for the constant kernel K = 2. This kernel is technically 
easiest to deal with, and the general framework is most transparent. The- 
orems 4.2, 4.3 and 4.4 are the main technical results and serve to establish 
Theorems 2.1 and 2.2 for this kernel. 



4.1 Preliminciries 

Smoluchowski's equation with constant kernel K = 2 has a unique global 

solution in an appropriate weak sense given any initial size-distribution mea- 
sure with finite zcro-th moment [16, §2]. For convenience, we adopt the 
normalization in (2.3). The moment identity (2.1) is valid for all bounded 
continuous functions f on E, and taking / = 1 we find that the total number 
density of clusters is utiE) = t~^. Since tvt{E) = 1, we associate to each 
solution a probability distribution function 

Ft{x)= [ i^tidx) / [ vt{dx) = tvt{x). (4.1) 

^(0,0;) / JE 

We also introduce the s-measures 

Gt{dx) = xvtidx), (4.2) 
and associated Laplace exponents 

<f{t,q)= [ {l-e-nMdx)= [ ^^^^Gtidx). (4.3) 

JE JE X 

Notice that q v{t,q) is strictly increasing with (p(t, 00) = ut{E) = 
t^^, and dq(p{t, q) is the Laplace transform of the mass-distribution measure 
X vt {dx) , so is completely monotone, if solves the simple equation 



(4.4) 
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for which the solution at any time t > is determined from data at time 
io > according to 

1^ = i^/^l°w, V q>0,t>0. (4.5) 

Since < (p(to,q) < t^^, we see that given Ff^ = toVto an arbitrary prob- 
abihty measure, (p{t,q) is well-defined on the time- interval (0,oo). But for 
< t < to, (p{t,q) may not have completely monotone derivative, and thus 
may not define a (positive) measure. The map q i-^ dqip{t, q) is completely 
monotone for all t G (0, oo) if and only if v is an eternal solution. 

Our study of convergence properties for solution sequences is based on 
pointwise convergence properties of 99, which are equivalent to convergence 
properties of the s-measures Gt according to the results of section 3. We 
begin by proving the continuous dependence of solutions on initial data, 
based on the evident fact that ip{t,q) is a continuous function of (p{to,q)- 

Theorem 4.1. (Continuous dependence on data.) For Smoluchowski's equa- 
tion with constant kernel K = 2, let to > and let i/^"^ be a sequence of 
solutions defined for t >to. 

(a) Ifi't!^^ converges weakly to a measure uq with i>Q{E) = t^^ , then for ev- 

(n) 

ery t > to we have that vl converges weakly to vt, the time-t solution 
with initial data vto = i'o- 

(h) For any t > to, if y^^^ converges weakly to a measure v with viE) = 
t~^, then v^^^ converges weakly to a measure uo with C'o{E) = t^^ , and 
u = ut, the time-t solution with initial data uto = i'o- 

Proof. We prove part (b); part (a) is similar. Let G^-j^\dx) = xu^'^\dx) 
and G{dx) = xi>{dx), and let (p^"'\t,q) and (p{q) be the associated Laplace 
exponents as in (4.3). The hypothesis is equivalent to saying that the s- 
measures G["'\dx) converge to a non-divergent s-measure G{dx) = xv{dx) 
with f^x~^G{dx) = t~^. This is equivalent to the statement that for all 
q > 0, f^"'Ht, q) — >• (f{q) as n ^ 00, where (p{oo) = t~^ and 0(0'^) = 0. Then 
it follows that 

Since (^o(O^) = and ipo{oo) = t^^, we conclude that (po is the Laplace 
exponent for a measure z>o on E with i>Q{E) = t^^ , and that u'f^'' converges 
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weakly to i>o. We compare (4.6) with the exphcit solution (4.5) to see that 
v = vt, where u is the solution on [to, oo) with initial data vq. □ 

4.2 The scaling attractor and eternal solutions 

We are ready to prove Theorem 2.1 for the kernel K = 2. First we con- 
sider part (b), the correspondence between the scaling attractor and eternal 
solutions. 

Theorem 4.2. A probability measure F is an element of the scaling attrac- 
tor for Smoluchowski 's equation with constant kernel K = 2 if and only if 
F = vi for some eternal solution v. 

Proof. Let us first suppose that F = vi for some eternal solution v and 
show that F G Ap. Pick arbitrary sequences Tn,l3n oo, and consider the 
sequence of rescaled eternal solutions 

i^t\dx) = ^i^t/TAl^n'dx), t>0. (4.7) 

Observe that i'^\e) = t^^, therefore, 

F^\Pnx) = u,{x) = F{x) 

for every x. Thus, F G Ap by Definition 2.3. 

Conversely, suppose F € Ap. Wc shall show that F = ui for some eternal 
solution V. Let (p denote the Laplace exponent of F, and L'^^\Tn, /3n be as 
in Definition 2.3. Consider the rescaled measures 

Pi")(ri.x)=T„4)(/3„dx). 

This rescaling yields a solution that is defined for t > 1/T„, and by hypoth- 

(n) 

esis we have that converges weakly to F. Then by Theorem 4.1, for any 

t > we infer that i/l^^ converges weakly to ut where vtiE) =t''^ and is a 
solution with ui = F. The solution u is eternal since it is defined for t>to 
for every to > 0. □ 

Let us now prove that Ap is invariant (part (a) of Theorem 2.1). Suppose 
is a solution on some time interval [ti,oo), normalized so t'j(-E') = t~^. 
Suppose Ft € Ap for some T >t\. Replacing vtidx) by TvTtidx), we may 
presume T = 1 without loss of generality. By Theorem 4.2 above, Ft = i>i 
for some eternal solution 9. But then vt = Vt for all t > ti, meaning that 
V is (the restriction of) an eternal solution. We obtain that Ft G Ap for all 
t > by a similar scaling argument. 
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4.3 Levy-Khintchine representation of eternal solutions 

Theorem 4.3. (a) Let v he an eternal solution to Smoluchowski's equa- 
tion with K = 2. Then there is divergent s-measure H such that as 
t [0, the mass measure Gt{dx) = xvfidx) converges to H. 

(h) Conversely, given any divergent s-measure H there is a unique eternal 
solution with the properties in part (a), defined for all t G (0, oo) via 

Proof. We first show (a). Immediately from the solution formula (4.5), 

lim ^{t, q) = hm ^ ^^^'^^\ , = =: $(,) (4.9) 

exists for all g > 0, with <^{q) < oo, $(0+) = 0, and <I>(oo) = oc. By 
Theorem 3.1, Gt converges to an s-measure H with Laplace exponent 
and is a divergent s-measure by the criteria in (3.5)-(3.6). 

Let us now prove (b). Let H he a divergent s-measure with Laplace 
exponent By (4.9), any eternal solution with the properties in part (a) 
must be determined by (4.8). Observe that the function q/(l + tq) has 
completely monotone derivative for t G (0, oo). It follows that dqip{t,q) is 
completely monotone when (p{t,q) is given by (4.8) [9, XIII.4]. Moreover, 
with Uf determined from (4.3), ft(£') = (p{t, oo) = t~^. Thus, Vf is indeed 
an eternal solution. □ 

Remark 4.1. Observe that Gt{E) = xvt{dx) is finite for some t G (0, oo) 
if and only if it is finite for all t. However, it is not necessary that the mass 
be finite for a solution to be well-defined. 

Theorem 4.3 establishes parts (a) and (b) of Theorem 2.2. To establish 
part (c), we need to show that the map vi H from Ap to is a bi- 
continuous bijection. 

Theorem 4.4. Let i^^"^ be a sequence of eternal solutions with corresponding 
divergent s-measures Fixt>{). Then, taking n oo, the following 

are equivalent: 

(i) Jf^^ converges weakly to some measure v with i>{E) = t~^. 
(ii) converges to some divergent s-measure H. 
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// either (equivalently both) of these conditions hold, then v = for an 
eternal solution with divergent s-measure H . 

Proof. Assume (i), so u\^^ converges to z> with i>(-E') = t"^. Then G^^\dx) = 
xv^^^dx) converges to G{dx) = xi>{dx) and the associated Laplace expo- 
nents converge: (p^"'\t,q) i^{q) for ah g > 0. Hence 

as n ^ oo for every q > 0. Since <^(0"'") = and tip{q) 1 as q oo, < 
oo for every q > 0, <&(0^) = 0, and hniq^oo ^(q) = oo- By Theorem 3.1 and 
(3.5)-(3.6), this proves (ii). 

We now show (ii) imphes (i). Suppose the divergent s- measures H'^"-^ 
converge to a divergent s-measure H. Then Theorem 3.1 with (3.5)-(3.6) 
imphes ^{q) for every q > 0, $(0+) = 0, and ^ oo as 

g ^ oo. Then, 

^'"'('•"'"TTM- iris 

(n) 

for every q > 0. This yields weak convergence of to ft, where v is the 
eternal solution with Laplace exponent <I> and divergent s-measure H. □ 

4.4 Scaling limits and initial tails 

We now prove Theorem 2.5 for the constant kernel. 

Proof of Theorem 2.5. Introduce rescaled solutions i'^'^\dx) = TnV^ {[3ndx), 

and let f/"^ = tDf'\ Also let Gf'\dx) = xi>i"'\dx) and let ip^''\t,q) be the 
associated Laplace exponent. Let H be the divergent s-measure correspond- 
ing to u and $ its Laplace exponent, and let ip{q) be the Laplace exponent 
of G{dx) = xi'i{dx). 

Then statement (i) of Theorem 2.5 is equivalent to saying — >■ F 
weakly, meaning the s-measures G^"^ converge to G with x~^G{dx) = 1. 
This is equivalent to saying 

(^(")(1, q) ip{q), q>0, where 99(0+) = 0, 99(00) = 1. (4.11) 

On the other hand, since Tnxv'^\l3ndx) = Gijx^{dx)., statement (ii) of 
Theorem 2.5 is equivalent to the assertion 

(^(")(T-\ q) $(g), g > 0, where $(0+) = 0, $(00) = oc. (4.12) 
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But by the solution formulae (4.5) and (4.8), we have 



l + (r„-l-l)^W(l,g)' I-Viq)' 
Since evidently (4.11) is equivalent to (4.12), (i) is equivalent to (ii). □ 

4.5 The representation at +oo 

For the additive kernel, Bertoin showed that an eternal solution can be 

uniquely identified by its asymptotic behavior as t ^ oo also. For the 
constant kernel, an analogous result follows easily from (4.3) and (4.8). 

Theorem 4.5. Let v he an eternal solution of Smoluchow ski's equation with 
constant kernel K = 2, and let $ he the Laplace exponent of the divergent 
s-measure associated with v. Then as t ^ oo, the measure t^vt converges 
weakly on (0, oo) to a measure A+ with Laplace transform 

^+(q) := / e-«^A+(da;) = — ^ = lim / Q-'^'^iytidx). 

Jo ^{q) Jo 

Clearly an eternal solution v is uniquely determined from A_|_ through 
^{q) = l/^+{q). We see that the measure A+ has a Laplace transform 
$+(g) defined for all q > 0, and <&+(g) ^ oo as g ^ since $(0+) = 0. So 
Jq A^{dx) < oo and J^A+{dx) = oo. 

The class of measures A-|- which arise in this way is characterized by the 
property that r]{q) = dq{l/^^{q)) is the Laplace transform of some divergent 
s-measure H (i.e., r] is completely monotone, locally integrable on [0, oo) and 
Je ^(^) ^'i ~ °°)- There does not appear to be a simple characterization by 
moment conditions. 

Remark 4.2. This representation has an interesting probabilistic interpre- 
tation; see [2, p. 74]. If X. is a subordinator with Laplace exponent then 
= 1/$ is the Laplace transform of the potential measure U, defined on 
Borel sets A c E hy U{A) = E {J^ l{x,eA}ds) . 



5 The additive kernel 

In this section we study the scaling dynamics for the additive kernel. Our 
main aims are to prove continuous dependence on initial data, establish the 
correspondence between points on the scaling attractor and eternal solutions, 
and revisit Bertoin's Levy-Khintchine representation with convergence of s- 
measures in mind. 
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5.1 Solution by Laplace transform 

The solution of Smoluchowski's equation with kernel K = x + y hy the 
Laplace transform is classical [8], and remains the basis for rigorous work. 
Let to £ be arbitrary. We assume f^,, is a (possibly infinite) measure with 
J^xftoidx) < oo. Without loss of generality, we may assume J^xvtoidx) = 
1. We have shown [16, Thm 2.8] that (1.1) has a unique solution i/f for t>tQ 
in an appropriate weak sense, such that 

/ xi^tidx) = 1, t> to- (5.1) 
Je 

As for the constant kernel, we use the notation 

ip{t, q)= [ {I- e-'i^)iyt{dx), q>0. (5.2) 
Je 

and set (Po{q) = ^ito,Q)- To study scaling limits we consider the mass 
distribution function, which is the natural probability distribution function 
associated to a solution. Let 

Ftix) = I yvt{dy). (5.3) 
Note that the Laplace transform of Ft is 

/ e-'i-Ft{dx)=dq^{t,q). (5.4) 
JE 

Thus, dqip{t,q) is completely monotone and dqip{t,0) = 1, t > 0. We know 
from [16] that if we substitute f{x) = 1 — e~'^^ in (2.1) we find that ^p{t,q) 
solves the hyperbolic equation 

dt(p - (fdqip = -(p. (5.5) 
Following Bertoin, it is convenient to introduce the new variables 

s = e\ so = e*", i;{s,q) = ^^-^(t,^y (5.6) 
By (5.1) and (5.2), ip is the Laplace exponent 

^{s,q) = I (e-« -l + qy) Giogs{dy), (5.7) 
Je 
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where Gt denotes the s-measure 

Gt{dx)=x'^iyt{e^dx). (5.8) 

Observe that Gt is not a finite measure in general, but if Gt{E) < oo for 
some t, then Gt{E) is finite for for every t for which the solution is defined, 
and is constant. We substitute (5.6) in (5.5) to see that V satisfies the 
inviscid Burgers equation 

+ i^M = 0' s > so. (5.9) 

The values of V', dqil) and are positive for s > sq; 9 > 0, and d'^il:{s^ ■) 
is completely monotone since it is the Laplace transform of Gt- In addition, 
(5.1), (5.7) and (5.8) imply 

Mm. dq'il){s,q) = I x''^G\o^s{dx) = . (5.10) 

We may describe ip{s,q) globally for s > sq by the method of charac- 
teristics. A surprising fact is that we may always solve for ^ backwards in 
time, for all s > 0, without developing singularities. The solution need not 
correspond to a positive measure for t < to, however. This is analogous 
to the situation for iiT = 2. 

Lemma 5.1. Let to R and ft,, G M with x utf^ldx) = 1, and let 
V'o('?o) = Qo/sq — ipo{qo/so). There is a unique solution ijj{s,q) to (5.9) 
defined for every s > and q > 0, such that ip{so, •) = ^o(')- 

Proof. Applying the method of characteristics as usual, the solution = 
^{s,q) is determined implicitly from the equation 

h{s, q, ip) :='(p- tpo{q - {s - so)V') = 0- (5-11) 

We have h{s,q,0) < 0, and d^h > s/sq since dqipo < (5.10). Since 

ipo is analytic, (5.11) determines a solution of (5.9) analytic in {s,q) for all 
s>0, q>0. □ 

Equation (5.11) determines the solution at time s from data at time sq 
and plays the same role in the analysis here as equation (4.5) played in the 
previous section. Convergence properties of solutions will be deduced from 
the pointwise convergence properties of the Laplace exponent tp using the 
theory from section 3. 
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Theorem 5.2. (Continuous dependence on data.) For Smoluchowski's equa- 
tion with additive kernel K = x + y, lettQ^M. and let u^"'^ be a sequence of 
solutions defined for t>to with J^xv^\dx) = 1 for all t > to. 

(a) If xvl^\dx) converges weakly to a measure xvQ^dx) with J^x{>o[dx) = 

1, then for every t > to we have that xv^\dx) converges weakly to 
xvtidx), the time-t solution with initial data = uo. 

(b) For any t > to, if xv^^\dx) converges weakly to a measure xv{dx) with 
J^xi'{dx) = 1, then xUf"\dx) converges weakly to a measure xi}o{dx) 
with J^xuoidx) = 1, and v = Vf, the time-t solution with initial data 
^0 = t'Q- 

Proof. We prove (a); the proof of (b) is similar. Let G^^\dx) = x'^ut{e^dx), 
and with s = e* let 

^(-)(s,g)= / y-\e-'iy-l + qy)Gf'\dy). (5.12) 

The family ^*^"^)(so, •) is uniformly Lipschitz, since equation (5.10) implies 
^H(so, 0) = 0, < SgV^^^so, q) < l/so for all q>0. (5.13) 

(n) 

The hypothesis is equivalent to saying that the s-measures G^^^ converge 
to a non-divergent s-measure Go{dx) = x^z>o((ix) with J^x^^Go{dx) = 1. 
By Theorem 3.1 and the criteria in (3.5)-(3.6), this is equivalent to the 

statement that for all q > 0, '(/;(")(so, g) i'oiQ): where ^/Jo is the (second- 
order) Laplace exponent for Go. with dg%ljci{0'^) = 0, 9g'(/'o(c>o) = I/sq. (Note 
dqip is the first-order Laplace exponent of Gq.) As in (5.11) we have 

^^"\s,q) - i^^^\so,q- (s - so)^^^\s,q)) = 0. (5.14) 

For fixed s, q, the sequence ip^'^^ (s, q) is bounded, and any subsequential limit 
must satisfy 

ip* - ipo{q - (s - so)ip^) = 0, (5.15) 

due to the equicontinuity of the maps ip t-^ ij/^^Sdjq — (s — so)V')- But 
equation (5.15) has the unique solution -0* = '^{s,q), where is the so- 
lution of (5.9) with ilj{so,q) = 'ipoil), q > 0. Hence the whole sequence 
ilj^'^\s,q) converges pointwise to ^j{s,q). Moreover, differentiating (5.15) 
yields 9gV'*(0) = 0, 9g^/'*(oo) = l/s, since sq = 1. Then the conclusion of 
(a) follows from Theorem 3.1, (3.5) and (3.6). □ 
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5.2 The scaUng attractor and eternal solutions 

Theorem 5.3. A probability measure F is an element of the scaling attrac- 
tor Ap for Smoluchowski 's equation with additive kernel K = x -\- y if and 
only if F{dx) = xvQ{dx) for some eternal solution v. 

Proof. Suppose F{dx) = xuQ{dx) for some eternal solution u. We show 
F E Ap. Pick arbitrary sequences Tn,bn — > oo, and consider the sequence 
of rescaled eternal solutions 

(dx) = b-'ut-Tr, {b-'dx) , teR. 

The corresponding distribution functions satisfy F^{bnx) = F{x) for every 

X. Thus, F G .4p by Definition 2.3. 

To prove the converse, suppose F G Ap. We show that F = Fq for some 
eternal solution u. Let ip correspond to u as in (5.2), and v^'^\Tn,bn be as 
in Definition 2.3. Consider the rescaled measures 

i'^'\dx) = bni'l'^j'n^bndx). 

This rescaling yields a solution that is defined for t > — r„. By assumption, 

^o^"^(^) = fythdy) = F^\bnx) - F{x), 
Jo 

at all points of continuity. By Theorem 5.2, this implies that for any AT G N 

(n) 

the solutions i^^ converge weakly to Vf for all t > —N. In particular, is 
a solution for t > — iV for all N, thus it is an eternal solution. □ 

Let us now prove that Ap is invariant (part (a) of Theorem 2.1). The 
proof is substantially the same as for K = 2. Suppose is a solution on 
some time interval [ti,oo), normalized so J^xi^tidx) = l,t > ti. Suppose 
Ft G v4p for some T > ti. We may presume T = without loss (if not, 
we translate in time, replacing i>t{dx) by vt-ridx)). By Theorem 5.3 above. 
Ft = xijQ for some eternal solution P. But then vt = vt for all t > ti, 
meaning that is (the restriction of) an eternal solution. We obtain that 
Ft G for every i G M by a similar argument. 

5.3 Levy-Khintchine representation of eternal solutions 

We now prove Bertoin's Levy-Khintchine representation for eternal solu- 
tions. The proof mainly follows [4], and is included to stress the basic 
framework. 
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Theorem 5.4 (cf. Bertoin [4]). (a) Letv he an eternal solution to Smolu- 
chowski's equation with K = x + y, and let Gt{dx) = x'^Vti'^^dx) he 
associated s-measures. Then there is a unique divergent s-measure H 
such that Gt converges to H as t —oo. 

(h) Conversely, given a divergent s-measure H there is a unique eternal 
solution with the properties in part (a), defined as follows. Let 

*(g)= f l^l^l±^H{dx) (5.16) 
Je ^ 

he the Laplace exponent of H, and let ip = ilj{s,q) he the solution to 

ip-^iq- #) = 0. (5.17) 
Then vt is determined by (5.7) and (5.8). 

Proof. We first prove (a). By Theorem 3.1 and (3.5)-(3.6), it is enough to 
show that ^'(g) := \lm.s^Qtp{s,q) exists for every <? > 0, with dq^{^) = 
and dq^{oo) = oo. We know ip > ^ and dqip > 0, so dsip{s,q) < for all 
q, s. Hence it suffices to show that for each q > 0, tp{s, q) stays bounded as 
s i 0. 

1. We first show ilj{s,q) stays bounded for q near 0. Choose gi > 
such that := qi — ■il){l,qi) = (p{0,qi) > 0. Then ip{s,q) = ijj{l,qi) along 
the characteristic line joining (0, g*) and (l,gi), so < ipi^-iQ) < V'(l)9i) 
whenever < s < 1 and < g < g*. (See Fig. 5.3.) 

2. For q > q* the complete monotonicity of g i-^- q~'^'ip{s,q) implies 
V'(s,g) < q^q^'^i){s,q*). 

3. We now show dq^{0) = and dq^{oo) = oo. Observe that solves 



^'(g) = V(i,g + *(g)), g>0. 



Therefore, 



1 -0,^(1, g + *(g))- ^^■'^> 

Since q) = q — ip{0, g), we have dqip{l, g) = 1 — dqip{0, g) — >■ as g — >■ 0, 
— > 1 as g ^ oo. Thus, 5q*(0) = and 5g*(oo) = oo. This proves (a). 

We now prove (b). Let H he a divergent s-measure and ^ be defined by 
(5.16). Note ^5^(0) = and ^^^'(oo) = oo by (3.5)-(3.6). Since ^^^(g) > 0, 
ip{s,q) is globally defined and analytic with ip{s,q) < q/s, and (5.9) holds 
for all s > 0, g > 0. With $ and ip defined by (5.6), (5.5) follows. 
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S 




Figure 5.1: Geometry of characteristics 



By the well-posedness theory in [16], we obtain an eternal solution through 
(5.4), provided we show that dqip{t, •) is completely monotone, which implies 
that it is the Laplace transform of a (positive) measure that we call xut{dx). 
Prom (5.17) we obtain that (p = 'p{t,q) satisfies 

q = ^ + ^{sip), (5.19) 

whence 

dg^ = ^7^^. (5.20) 

l + s*'(s(y9) ^ ^ 

Since q i— > 1 + s^'{sq) is positive with completely monotone derivative, the 
map g 1-^ (1 + s"^' (sq))^^ is completely monotone [9, XIII. 4]. We then infer 
that dq<f{s, •) is completely monotone by Lemma 5.5 below. Since ^'(0) = 
we have the normalization (5.1), dqip{t,0) = 1. This finishes the proof of 
existence. 

Note that total number of clusters vtiE) = (p{t, oo) = oo always here. 

Let us show that the eternal solution defined by this procedure is unique. 
Let i7 be a divergent s-measure and suppose u, u are two eternal solutions 
with s-measures Gt, Gt that converge to H. But this is equivalent to point- 
wise convergence of 'il){s,q) and ip{s,q) to ^{q) as s ^ where ip and 
solve (5.9). But the solutions to the inviscid Burgers equation with increas- 
ing initial data are unique, thus ^{s, q) = tp{s, q) and v = v. □ 
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Lemma 5.5. Suppose f,g: E ^ E, f = g{f) and g is completely mono- 
tone. Then f is completely monotone. 

Proof. We prove by induction that the first n derivatives of G o / alternate 

in sign for every completely monotone function G. For n = 0, G{f) > 
0. Suppose the statement is true for some ra > 0. Let G be completely 
monotone, and note 

-{Gofy = -G'{f)g{f) = G{f) 

and G is completely monotone since it is the product of completely monotone 
functions. Using the induction hypothesis, we deduce that the first n + 1 
derivatives of G o / alternate in sign. □ 

To complete the proof of Theorem 2.2 for K = x + y, we need to check 
that the map H from Ap to Sd is a bi-continuous bijection. 

Theorem 5.6. Let i^*-") be a sequence of eternal solutions with corresponding 
divergent s-measures Fixte^. Then, taking n oo, the following 

are equivalent: 

(i) xvf^^ converges weakly to xu with J^xv{dx) = 1. 

(ii) converges to a divergent s-measure H. 

If either (equivalently both) of these conditions hold, then v = vt for an 
eternal solution with s-measure H. 

Proof. With Theorem 3.1 in hand, the proof of Theorem 5.6 is essentially 
the same as that of Theorem 5.2. Assume (i), so u^"^ converges to u 
with J^xi>{E) = 1. Then G["\dx) = x'^v\^\e*dx) converges to the s- 

mcasurc G{dx) = x'^u^e^dx) and the associated Laplace exponents converge: 
i^{q) for all q > 0, with ^^^(O) = 0, dgij^ioo) = 1/s. Recall that 
tp^'^\s,q) solves 

= ^W(s^g). (5.21) 

Let M > 0. A calculation as in (5.18) shows that dq'^'^^''\q) is uniformly 
bounded in n for g G [0, M]. We claim that lim„_KX) ^^"^(g — s^(q')) exists 
for every q. Let us restrict attention to g G [0, M]. Then by (5.21) 

^(n) (^ _ g^(^)) ^ ^(n) ^) ^ j^^{n) _ ^^^^y^ _ ^(n) _ ^^(n) ^^^^ ^ 

The first term converges to ip{q) and the second to zero by the uniform 
estimate on aq*(")(g) on [0,M]. Since M > was arbitrary, we may use 



G. Menon and R. L. Pego 



33 



Theorem 3.1 to deduce that '^^'^\q) converges to a Laplace exponent ^{q) 
that satisfies 

^{q-si^{q)) =i>{q). 

As with (5.18) and its sequel it follows that dq"^{0) = and dq'i^{oo) = 
oo. Thus ^ is the Laplace exponent of a divergent s-measure H, and 
converges to H. 

We now show (ii) implies (i). Suppose the divergent s- measures H^""^ con- 
verge to a divergent s-measure H. Then Theorem 3.1 implies ^("^(g) ^(q) 
for every q > 0, and (9q^(0) = 0, dq'^{oo) = oo. Then the characteristics 
emanating from s = converge because q + s'^^'^\q) q + s^{q). Thus, 
tp^'^\s,q) tp{s,q), which satisfies (5.17). This yields weak convergence of 

(n) 

xvl to xvt, where v is the eternal solution with divergent s-measure H. □ 

5.4 Scaling limits and initial tails 

Let us now prove Theorem 2.5 for the additive kernel. 

Proof of Theorem 2.5. We rescale solutions via i'\^\dx) = (5nV^j,^{l3ndx) , 

and let Fi''\dx) = /3„xz>^"^((ix). Also let G^^\dx) = x'^v[''\c*dx) and let 
ijj^^\s,q) be the associated Laplace exponent as in (5.7). Observe G^") m 
(2.24) is <5^'^^ and J^x~^&"'\dx) = e^". Let H be the divergent s-measure 
corresponding to u and * its Laplace exponent, and let ip{q) be the Laplace 
exponent of G{dx) = x^uo{dx). 

~ in] 

Then (i) is equivalent to saying Fq ^ F weakly, meaning the s- 
measures G'q'^ converge to G with J^x~^G{dx) = 1. This is equivalent 

to saying 

^(")(1, q) ip{q), q>0, where dg'ip{0+) = 0, dgtpioo) = 1. (5.22) 

On the other hand, since /3„x^fQ"'^(e~-^"/3„dx) = G-Tn{dx), (ii) is equivalent 
to saying 

?/;(")(e~^",g) ^ *(9), q>(}, where aq*(0+) = 0, 5g*(oo) = oo. 

(5.23) 

For brevity, let ip^'^'^\q) denote V'(")(l,g) and ^^^'^q) denote ■i/'(")(e~^", g). 
Then the implicit solution formulas to (5.9) read 

^(g) = *(g - V'(g)), = *^"^(« - (1 - e-^")V'^"^(g)). 

As in the proof of Theorem 5.6 we may now deduce that (5.22) is equivalent 
to (5.23), implying (i) is equivalent to (ii). The details are omitted. □ 
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6 The multiplicative kernel 

In this section we study scaling dynamics approaching the gelation time for 
the kernel K = xy. The study of the multiplicative kernel can be reduced to 
the additive kernel by a simple change of variables. This trick is well-known 
(see [8]), and allows us to avoid separate proofs. 

6.1 Solution by the Laplace transform 

The self-similar solutions for K = xy have infinite number and mass, but 
finite second moment. However, one may develop a natural well-posedness 
theory using only the finiteness of the second moment [16]. We assume 
is a (possibly infinite) measure with x'^uto{dx) < oo. Without loss of 
generality, we may scale so that J^x^vtoidx) = 1 and to = —1 as in (2-3). 
We define the Laplace exponent (note the change from (5.2)) 



Je 

and write ipo{q) = (p{to,q). We may substitute (6.1) in the moment identity 
(2.1) to obtain 



Equation (6.2) may be transformed to (5.5) by the following change of vari- 
ables. Let ip^'^'^{T,q), T S [0,oo), denote a solution to (5.5) with initial data 
</?o(<z)- Then the solution to (6.2) is given by 




(6.1) 



dtip - ipdqip = 0, te {to, 0). 



(6.2) 



^(t,^)=eV^''(r,g), r = log(|t|-i), 



(6.3) 



which may also be written in terms of the number density as 

xutidx) = e'^uf'^idx). 



(6.4) 



Conservation of mass (5.1) is now replaced by 




(6.5) 



and the probability measure Ft associated to vt is defined by 




(6.6) 



As in (5.8) we define the s-measure 

Gt{dx) = x^uti\t\-^ dx) = Gf'^idx), t G [to, 0), 



(6.7) 
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and the associated Laplace exponent 

i^{t,q)= I y-^{e-'iy-l + qy)Gt{dy)=i;^''{\t\-\q). (6.8) 
JE 

The correspondences (6.4), (6.7) and (6.8) map normahzed solutions for K = 
xy on the time interval t G [—1,0) to normalized solutions with K = x + y 
on the interval r G [0, oo). The same change of variables may be applied 
to eternal solutions defined on t G {—oo,0). By consequence, the results 
established so far for the additive kernel carry over in an obvious way for 
the multiplicative kernel. This yields continuous dependence of solutions on 
data (by Theorem 5.2), the correspondence between the scaling attractor 
and eternal solutions (Theorem 2.1), the Levy-Khintchine representation 
(Theorem 2.2), and how initial tails encode scaling limits (Theorem 2.5). 
For completeness, we make explicit the map from divergent s- measures to 
eternal solutions implicit in Theorem 2.2(b). 

Theorem 6.1. Given a divergent s-measure G there is a unique eternal 
solution defined as follows. Let 

m = I ^^^^^^±^H{dx), 

Je 

and 'tp{t,q), t £ (— oo,0) be the solution to 

ip - ^'(g + t"V) = 0. 
Then is determined by (6.7) and (6.8). 

7 Doeblin solutions 

This section is inspired by Feller's treatment of Docblin's universal laws 
and domains of partial attraction [9, XVII.9]. But apparently we must be 
content with using more words to prove fewer results. Our aim is to prove: 

Theorem 7.1. There exists an eternal solution v whose scaling ui-limit set 
contains every element of the proper scaling attractor, . 

We will show later that A is the closure of (see Corollary 9.5). There- 
fore, Theorem 7.1 establishes Theorem 2.6. 

The proof is based on suitably "packing the tails" of the corresponding 
divergent s-measure. The following is adapted from Feller [9, XVII.9]. Given 
an s-measure G and a, 6 > we define a rescaled measure G"'^ by 

G'''\x) = aG{bx). (7.1) 



(6.9) 
(6.10) 
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Lemma 7.2. Let he a sequence of s-measures with 

x-^Gk{dx) < k. (7.2) 



Ie 



oo 



IE 

Then there exist sequences ak,bk such that ak 0, a^bk oo, 

k=l 

defines an s-measure, and — Gk converges to zero. 



G:=^Gl'"''' (7.3) 



The growth assumption (7.2) is included only for concreteness and im- 
plies no real loss of generality. Our main purpose is to approximate divergent 
s-measures. 

Lemma 7.3. Let H be a divergent s-measure. There exists a sequence of 
s-measures G^ satisfying (7.2) such that G^ converges to H . 

Proof of Theorem 2.6. 1. Let u^'^^ be an arbitrary sequence of eternal solu- 
tions with corresponding divergent s-measures H'^'^\ Partition the integers 
into infinitely many subsequences, and choose Gk H^'^^ for k in the n-th 
subsequence as in Lemma 7.3. 

2. Now define a^, h^ and G as in Lemma 7.2, and put 

H = 5q + G. 

H has an atom at the origin, thus is the divergent s-measure for an eternal 
solution. By construction, G"''"^'' H^""^ as k ^ oo along the n-th subse- 
quence. Moreover, since — > 0, under rescaling S^'"^'' = UkSo converges to 
zero. Thus, if we take limits along the n-th subsequence, H°'k,bh jjin) 

3. We now apply Theorem 2.2 together with (8.1). We have = 
6p(F,";''"=) and (x) = FT,{f5kx) where 



{Tk,(ik) 



{akbk,bk) {K = 2), 

{\ogiakbk),akbl) {K = x + y), (7.4) 
, {-{akbk)~^,akbk) {K = xy). 



Observe that Tk —> Tmin and 0k 'Oo. We take limits along the n-th 
subsequence to obtain FT^[f3kx) F^''^\x) at every point of continuity. 
Hence, given any sequence of eternal solutions v^"'^ there exists an eternal 
solution V whose scaling co-limit set contains each I'l ' . 
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4- The space of divergent s-measures is separable. The s-measures which 
are concentrated at finitely many rational points (including 0) with rational 
weights form a countable set which is dense with respect to convergence of 
s-measures. By ordering these in a sequence ^(") and using the construction 
above, we see that there exist eternal solutions u such that for every eternal 
solution u, ui is in the scaling a;-limit set of u. This finishes the proof of the 
theorem. □ 



7.1 The packing lemma 

We will need to choose a sequence Ck that grows so fast that 

oo 
j=k+l 

1,2 

The choice Ck = e will do. For j >2 we have the elementary estimate 
je~^ < / ye~'^ dy = cosh j. 

Jj-l/2 

Therefore for > 1, 



2 A .2 .2 „2. e-'^-V^ 



J2 Je ^ < 6^= / ^ ye ^ dy = > 0. 

j=k+i 



k+l/2 2 



Proof of Lemma 7.2. 1. Fix akbk = c^. Then G defines an s-measure since 

» oo „ oo 

/_ x-'^G{dx) = ^c^^ /_ x-^Gkidx) < ^ kc^^ < oo. 
•^^ k=i •'^ k=i 

2. Let ^^'^^ and $ denote the Laplace exponents of G^ and G respectively. 
We use the definition (7.1) and (7.3) to obtain $(g) = Ejli cj^^^-^'^^^'j)- 
Observe that G"*" — Gjt is a positive measure with Laplace exponent 

ck^qb-') - $«(g) = Cfe5^cTi$0)(g6,.6-i). 

3. To prove convergence to zero, it suffices to show that the right hand 
side converges to zero for every q > 0. We first control the tail. Since 
JgX-^Gj{dx) <j, 

oo oo 
j=k+l j=k+l 
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4- We now choose bk inductively to control the first A; — 1 terms in the 
range < q < k. Suppose bi, . . . , b^^i have been chosen. Since 
as g ^ 0, we choose b^ so large that 

1 '"'^ 1 
ak = Ckb,' < Cfe^c7^$(^) {kbjb^') < -. 

□ 



7.2 Proof of Lemma 7.3 

First, suppose H has no atom at the origin. Since y^^H{dy) — > as 
X — oo, we may choose a decreasing sequence such that y~^H{dy) < 
k. Let Gk{dy) = ly>ei,H{dy). Clearly, Gk satisfies both conditions of Defi- 
nition 2.5. 

Next, let H = Sq. In this case we choose Gk{dx) = {xlogk)^^l^>k-idx. 
Then Gj- satisfies (7.2) as 

/ x-^Gkidx) = (logk)-^ / x-^dx = k{logk)-^ < k, 

and 

= q(logky^ / — dx q = ^{q). 

The general case follows by superposition of these two special cases. 



8 Scaling-periodic solutions 

In this section we characterize scaling-periodic solutions and show that they 
are dense in the scaling attr actor. That is, we prove Theorems 2.7 and 2.8. 



8.1 Characterization 

Proof of Theorem 2.7. 1. Given a scaling-periodic solution, a solution sat- 
isfying (2.25), we can scale it as in (2.16)-(2.17) so that to is as in (2.6). 
Then, under the map F i— ^ F"''^ given by 



Fabtibx) {K = 2), 

Ft+iog{ab) {ab'^x) {K = x + y), 
Ft/ah{ab^x) {K = xy), 



(8.1) 
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for some a, 6 > we have Ft = F"'^ for all t G [to, T^ax)- Explicitly, 
{ab, b) 



' {h,P) {K = 2), 

(e*i,/3e-*i) {K = x + y), (8.2) 
^ {K = xy). 



Observe that a6 > 1 in all three cases. Iterating the map, we get that the 
solution must be (the restriction of) an eternal solution. By Theorem 2.3 
and (8.1), F = F"'^ is equivalent io H = that is, 

H{x) = aH{bx), x > 0. (8.3) 

Without loss of generality we may suppose 6 > 1 since (8.3) is equivalent to 
a-^H{b-^x) = H{x). 

2. Equation (8.3) implies -ff(0+) = ai/(0+). If H has an atom at the 
origin, this forces a = 1. Then H{x) = H{bx) for every a:; > 0, and since 
6 > 1 and H{x) is non-decreasing, it follows H{x) = c = i?(0+) for all x > 0. 
Therefore, if H has an atom at the origin, then H = cSq for some c > 0. 

3. Suppose H does not have an atom at the origin. We iterate (8.3) to 
find that 

H(dx) , f-' H(dx) 



In order that H is an s-measure we require 
/ (1 A x-^)H{dx) = J2 / ^id^) + J2^ab)-^ j x-^H{dx) < DO. 

Thus, a < 1 and ab > 1. Given a; > let A; = max{j : V < x}. A similar 
calculation yields 

H{x) = Va-^' / H{dx) + a-'' [ H{dx). 
j<k 

This shows H is determined by its restriction to [1,6). 

4- Conversely, suppose H = H""'^ and (i) or (ii) hold. Notice that H is 
automatically divergent since it either has an atom at the origin or 

/ x-^H{dx) = / x-^H{dx) V {ab)-^ = oo. 

j=-oo 

Thus, it determines an eternal solution, which by (8.1) satisfies F = F""'^ . □ 
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8.2 Self-similar solutions 

As remarked in Section 2.6, the case (i) is simple but important. The asso- 
ciated divergent s-measure is scale-invariant for every b> 1 and the scaling- 
periodic solutions are the classical self-similar solutions with exponential 
tails. If a scaling-periodic solution satisfies (2.25) for every ti > to (with 
changing /?), it follows that for some fixed a and b, H{x) = aJ'H{Wx) for all 
rational and hence all real r. The fundamental rigidity lemma for scaling 
limits [9, VIII. 8] then implies H{x) = Cgx^ for some G M. The finiteness 
condition /g(l A x~^)H{dx) < oo then implies 9 = 1 — p, p & {0,1]. If p = 1, 
H is an atom at the origin corresponding to (i) above. The self-similar pro- 
files and their domains of attraction are discussed further in Section 10.2. 

8.3 Density of scaling-periodic solutions 

To prove Theorem 2.8 and establish density of scaling-periodic solutions in 
the full scaling attractor A, it will suffice to prove such solutions are dense 
in the proper scaling attractor Ap (see Corollary 9.5). 

Theorem 8.1. Scaling-periodic solutions are dense in Ap. 

Proof. 1. Let F ^ Ahc arbitrary. Let a„ J, 0, 6,i j oo be sequences such 

1/2 

that Oribn and ffln^n oo. We claim that there exist scaling-periodic 
solutions J^^"^ with scale parameters (a„, 6„) such that Fj^"'(d.T) = x'^v^1^\dx) 
converges weakly to F as n ^ oo. Let H denote the divergent s-measure 
associated with u. By Theorems 2.2 and 2.3 it suffices to construct divergent 
s-measures iJ^") such that H^""^ = anH^'^\bn ■) and //("^ converges to H. 

2. Consider first the case where H has no atom at the origin. In this case 
we define i?*^"-* to be the scaling-invariant extension of H restricted to the 
interval := ^n^^'^ ,b]i'^). Then for any a; > that is a point of continuity 
of H, ioT n large we have x E In and 




as n 



oo. Moreover, 



/' 

Jx 



oo 



Hi'^Xdy) 




'OO 



H{dy) 



y 



y 



This establishes the desired convergence of s-measures. 
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3. In case H = 60, we let H^'^^ be a sum of delta masses Sj^\j G Z 
concentrated at points = b'n'^^^, so that H^^^ = X]j(a„6„)-^5j"\ Observe 
that there is no mass in {bn^^^ ,bn'^); thus for any x > 0, for n large we have 

H(^){x) = Y^ianbuY = , ^ , ^ 1, 
j<0 

and 

Hence the s-measures i?^") converge to b^. 

4. In the general case, we simply superpose the separate constructions. 

1/2 . 

Observe that the restriction a„6„ — > is only needed in the critical case 
when H has an atom at the origin. □ 



9 Extended solutions, with dust and gel 
9.1 Extended solutions 

A proper solution to Smoluchowski's equation satisfies J^x'^Vfidx) = m^{t) 
with m^{t) normalized as in (2.4). However, a sequence of proper solutions 
may lose mass in the limit. We append atoms at and 00 to account for 
these defects, considering measures on = [0, 00] of the form 

Ht = ao{t)5o + x'^vt + aoo{t)5oo, (9.1) 

where vt is a size-distribution measure on E. We call the atoms ao and 
the dust and gel respectively. An associated probability measure on E is 
defined as in (2.5), by 

The s-measure associated to a solution in (2.14) is replaced by the s-measure 

Gtidx) = x^+'niXit) dx) + gooit)5^{dx), g^t) = (9-3) 

The measures /xt define Laplace exponents by evident modification of equa- 
tion (4.3) for K = 2, namely 

ifit, q)= [ {I- e-'i^)^it{dx) = aoo(t) + / (1 - e-^^^dx), (9.4) 
Je Je 
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and of (5.2) and (5.7) ioi K = x + y and (6.1) for K = xy, both yielding 
fit, q) = / fitidx) = ao{t)q + / (1 - e-«^)i/t(dx). (9.5) 

The evolution equations for these exponents remain 

-^^ (i^ = 2), 
dtip=l ifdgip-ip, {K = x + y), (9.6) 
ipdqif, {K = xy). 

This motivates the following definition. 

Definition 9.1. A family of triples (fj, ao(i); aoc(i))' ^ ^ [^Oj^max), defines 
an extended solution to Smoluchowski's equation for the kernels K = 2, x+y 
and xy with initial data {p, cq, cioo), if 

(a) The measures in (9-1) satisfy /xt(-E) = m^{t) with m^{t) as in (2.4), 

for t € [to,Tmax)- 

(b) (9.6) holds for g > and t G {to,T^^). 

(c) fxt ^ jl = aoSo + x'^v + 000(^00 weakly as f | to- 
Due to the normalization in (a), we regard extended solutions as deter- 
mined by the associated probability distributions F in (9.2). We will usually 
denote an extended solution with values (t'f, ao(t), aoo(t)) simply by v. 

Extended solutions provide the correct compactification in light of the 
following theorem. Since every proper solution also defines an extended 
solution, the theorem applies in particular to sequences of proper solutions. 

Theorem 9.1. Let Ff^\ t G [ioj^max); be probability measures associated 
with a sequence of extended solutions Then there exists a sequence 

rij 00 and probability measures Ft associated with an extended solution v, 

such that F^^^^ converges weakly to Ft for every t G [ioj^max)- 

Proof. Consider the sequence of probability measures F^^ ' on. E. Then 

there exists a subsequence nj and a probability measure Fq such that F^^ ^ 

converges weakly to Fq. We use Fq to determine initial data to define an 
extended solution z> for t G [^Oi^^max)- Continuous dependence on initial 

-(n) 

data as in Theorem 9.2 below implies the weak convergence of F^ ^ to Ft 
for every t G [to,7'max)- □ 
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We state the following result without proof, as it is an easy consequence 
of Definition 9.1, and Theorems 4.1 and 5.2. The notion of extended solution 

allows us to simplify matters, as it is no longer necessary to assume that 
fi{}{E) = m^(to) , or /i(-E) = m^(t) as in parts (a) and (b) of Theorem 4.1 
and 5.2. 

Theorem 9.2. ( Continuous dependence on data.) For Smoluchowski's equa- 
tion with kernels K = 2, x + y or xy, let to G (Tmin, ^max) o,nd let F^"^^ 
determine a sequence of extended solutions defined for t E I = [to, Tmax)- 

(a) If Ff^^ converges weakly to a measure Fq, then for every t E I, F^"'^ 
converges weakly to Ft, associated with the time-t extended solution 
with initial data determined by Ft^ = Fq. 

(b ) For any t E I, if F^^^ converges weakly to a measure F, then F^^^ 
converges weakly to a probability measure Fq and F = Ft, associated 
with the time-t solution with initial data determined by Ft^ = Fq. 

9.2 Transformation to proper solutions 

Clusters of "zero" or "infinite" size interact with other clusters in simple 
ways. The invariances of the evolution equations (9.6) allow us to relate all 
extended solutions (except pure dust and gel) to proper solutions. Let us 
consider the constant and additive kernels in turn. 



9.2.1 The constant kernel 

The dust and gel are recovered as limits as q ^ and cx3 respectively: 

aoo(i) = 'P{t, 0+), ao(i) = Ht{E) - ^{t, oo'). (9.7) 
Since Ht{E) = t~^, we take limits in (9.6) to see that the dust and gel satisfy 

The extended solution corresponds to purely dust and gel when utiE) = 0, 
so that ao{t) + aoo{t) = t~^. We may exploit (9.6) to show that every 
extended solution that is not purely dust and gel is in correspondence with 
a proper solution after a simple change of scale. Suppose ip{t, q) is the 
Laplace exponent of an extended solution. If aoo(io) > let 

0, q) = a{t)-^ {^{t, q) - a^{t)) , (9.9) 
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where 

£-1 = aitr\t-' - ao{t) - a^{t)), a{t) = (9.10) 

Then we find (p{i, 0+) = 0, (p{i, oo^) = i^^, and d^(p = —(p^. Thus, (p{i, q) is 
the Laplace exponent of a proper solution defined on [to, oo). 

For vanishing gel (aoo(to) ~^ 0) the transformation above simplifies, 
yielding a = l,t—t{i = t—tQ, = ip. Zero-size clusters combine trivially with 
other clusters, so the presence of dust only shifts time in accord with our 
normalization of total number. Observe that if gel is present (aoo(io) > 0), 
the probability of being gel approaches one {aoc{t) / fJ-tiE) 1) and the rel- 
ative distribution of finite-size clusters approaches a state reached by the 
proper solution at a finite time; we have t ^ to + l/aoo(io) as i — > oo. 



9.2.2 The additive kernel 

In this case, dq(p{t,q) = e""^^ ^t{dx) and iJ-t{E) = 1, so the dust and gel 
are given by 

ao{t) = dgifit, oo), a^{t) = 1 - dgifit, 0+). (9.11) 

The similarity with the constant kernel is clear if we use the time scale s = e* 
and the Laplace exponent ijj{s,q) defined in (5.6) and (5.7). Let 

6o(s) = ^-5,V(s,oc) = ^, b^{s) = dg^{s,0) = ^^. (9.12) 
We then take limits in (5.11) to see that 

d(s-l-6o) 1 2 ^^oo ,2 /niQ\ 

-^s = -^' ^ = ^^-'^^ 

which is equivalent to the following closed equations for the dust and gel: 

^ = -oo(l-ao), = aoo(l - floo)- (9-14) 

The extended solution is purely dust and gel when ao{t) + a^oit) = 1. If it 
is not, we exploit the invariances of the inviscid Burgers equation (5.9) to 
reduce extended solutions to proper solutions by a change of scale. Given 
initial data V'o with dq^po{0) = boo{so) > and dqipo{oo) = Sq ^ — 6o(so) > 
boo{so) > 0, we define a proper solution via the change of variables 

^P{s, q) = a{s)-^ {^{s, q) - b^{s)q) , (9.15) 
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where 



s ^ = a{s) ^ (s ^ - bo{s) - boois)) , q = a{s)q, a{s) 



boo{so)' 

This ensures dq'ip{s, 0) = 0, dqip{s, oo) = s~^, and dgip+i^dqip = for s > sq. 



9.3 Levy-Khintchine representation 

Definition 9.2. An extended solution to Smoluchowski's equation that is 
defined for all t G (Tmin, ^max) is called an eternal extended solution. 

The following representation theorem is the completion of Theorem 2.2. 
We establish a bijection between the set of eternal extended solutions and 
the space S consisting of all s-measures together with a point at infinity. The 
point at infinity corresponds to all measures such that f A x~^)H{dx) = 
oo. These measures give rise to the (unique) Laplace exponents ^{q) = 
'^{q) = oc, q > 0. We say a sequence of s-measures converges to the point 
at infinity if oo, q > for the associated Laplace exponents. This 

special case corresponds to the eternal extended solution that is pure gel. It 
is the counterpoint to the Laplace exponents $(g) = '^{q) = 0, q > which 
generate the eternal extended solution that is pure dust. 

Theorem 9.3. (a) Let u he an eternal extended solution for Smoluchowski's 
equation with K = 2,x + y or xy. If v is not pure gel, there is an s- 
measure H such that Gt converges to H as t I Tmin- If i' is pure gel, 
Gt converges to the point at infinity in S. 

(b) Conversely, for every s-measure H, there is a unique eternal extended 
solution u such that Gt converges to H as t I T^^a- The point at 
infinity generates the extended solution v that is pure gel. 

(c) Let & : A ^ S map the (full) scaling attractor A to the set S of s- 
measures by 6{F) = H, where H is the 's-measure associated to the 
eternal extended solution v such that F = Ft^ with to as in (2.3). 
Then & is a bi-continuous bijection. Moreover, &p: Ap — ^ is the 
restriction of & to Ap. 

The map & is defined in terms of Laplace exponents by the same formulas 
as for proper solutions: (4.8) for K = 2, (5.11) for K = x + y, and (6.10) 
for K = xy. Parts (a) and (b) of the theorem are then proven just as in 
Theorems 4.3 and 5.4. The proof here is simpler, since we no longer need 
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verify the divergence conditions on the s-measure. The proof of part (c) 
rehes on two separate arguments. It is easy to show as in Theorems 4.4 and 
5.6 that the map H i-^ Ftg is a bi-continuous bijcction. However, we must 
also identify such Ft^ as belonging to the attractor. Here the arguments 
deviate slightly from those of Section 4 and 5. We use parts (a) and (b) of 
Theorem 9.3 in the proof of part (c) via the following intermediate theorem. 

Theorem 9.4. (a) The scaling attractor A is invariant: If v is an ex- 
tended solution of Smoluchowski 's equation, and Ft E A for some t, 
then V is eternal and Ft ^ A for all t G (Tmin, Tmax)- 

(b) A probability measure F on E belongs to A if and only if F = Ft^ for 
some extended eternal solution v. 

Proof. The proof differs from earlier arguments only in the first part of The- 
orems 4.2 and 5.3 (the assertion that Ft^ E ^ if is eternal). In order to 
prove this, let us suppose v is an extended eternal solution with associated 
s-measure H = {H,g^) where H is an s-measure and g^o is the charge of 
y~^H{dy) at oo. To show that F := Ft^ is in the scaling attractor, we 
must find T„ | Tmax, f3n ^ oo and a sequence of proper solutions such that 
F!p\f3nx) F{x) at points of continuity. We use the Levy-Khintchine for- 
mula to find such solutions. We approximate H by the sequence of divergent 
s-measures H^"^^ = n~^do + H+goonSn- It follows that for the corresponding 
(proper) eternal solutions z/^"), the probability measures F^^^ converge to 
Ft for every t > T„^i^. Given any sequence T„ | T^i^i, Pn — > oo we consider 
a sequence of rescaled solutions determined as in (2.17), by 



We then have Fj^ {/3nx) = F^^ [x) —>■ F{x) at all points of continuity. 

The converse implication and part (a) are proven exactly as in Theo- 



This also proves a property alluded to several times before. 
Corollary 9.5. A is the closure of Ap. 

Proof. If F £ A has s-measure H, we approximate ^ by a sequence of 




Fljl{(3-^x), {K = 2), 
F^Xip-^x), {K = x + y) 
^t/lkl(/^n'(^)), {K = xy). 



f/")(x) = 



rems 4.2 and 5.3 and the sequel. 



□ 



divergent s-measures as above. 



□ 
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9.4 Scaling limits and initial tails 

We now state the natural extension of Theorem 2.5 to eternal extended 
solutions. The proof is almost identical to that of Theorem 2.5 except that 
we no longer need verify divergence of the s-measure. 

Theorem 9.6. Let F E A with associated s-measure H. Let be any 
sequence of proper solutions defined for t G [to^Tum^i) , with associated initial 
s-m,easures given by G'^^\dx) = x'^^^ v[^\dx) . Let Tn —<■ T^aax, Pn ^ oo. 
Then the following are equivalent: 

(i) Fj^\Pnx) F{x) as n ^ oo, at every point of continuity. 

(a) The resettled initial s-measures defined by (2.24) converge to the 
l-measure H as n ^ oo. 

10 Initial tails and ultimate scaling dynamics 

In this section, we present two applications of the principle that ultimate 

scaling dynamics arc encoded in the initial tails (as formalized in theo- 
rems 2.5 and 9.6). The first is a proof of the shadowing theorem 2.9. The 
second is a streamlined proof of the classification of domains of attraction 
in [16] that avoids the use of Karamata's Tauberian theorem. 

10.1 Initial tails and shadowing 

Proof of Theorem 2.9. 1. As in section 2, we let dist(-, •) denote any metric 
on V which induces the weak topology. Suppose that for Smoluchowski's 
equation with kernel K = 2, x + y or: xy, v and v are two solutions defined 
on [ioj^max), and make the assumptions stated in the theorem. Suppose 
that (2.30) fails, i.e., that 

dist(Ft(6(t) dx), Ft(bit) dx)) as t ^ r^ax- (10.1) 

Then since V is compact, by passing to subsequences we can find sequences 
Tn T Tinax and Pn = b{Tn) and different probability measures F, ]^ E P, 
such that as n — > oo we have 

FtAM ^ F{x), FfJPnx) ^ F{x), (10.2) 

at every point of continuity of the limit. Here the values T„, /?„ are those 
that correspond via the map {t, b) i— > (t, h) stated in the theorem. Relabeling 
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if necessary, we may assume < F{x) for some finite x, i.e., F does not 
represent pure gel. Therefore, according to the extended Levy-Khintchine 
representation theorem 9.3, there exists an s-measure H that corresponds 
to F. 



2. Let 



ar. 



f3n iK = 2), (Tn {K = 2). 

Pne-^" {K = x + y), Xn=\e^- {K = x + y), (10.3) 



Jn\Tn\ iK = xy), {\Tn\ {K = xy) 



and similarly define a^, in terms of /3n, Tn. Note that a„ — an- We 
claim that a„ — > oo. This is evident for = 2, and once we prove it for 
K = X + y it follows for K = xy hy the transformation formula (6.3). For 
K = x + y, one can prove q„ — ^ oc by following the beginning of the proof of 
Theorem 7.1 in [16] up to (7.8) using only subscqucntial convergence. Prom 
(7.8) one deduces Ae~* — >■ oo, which corresponds here to a„ oo. 
3. Define rescaled initial s-measures (sec (2.24)) by 

G(")(a;) = \nOi-^G{anx), &''\x) = \na-^G{anx), (10.4) 

According to the extended encoding theorem 9.6 the s-measures con- 
verge to H. Let f^^"^, (p^'^'^ and <I> be the first-order Laplace exponents associ- 
ated to G^'^\ G^^^ and H respectively as in (3.1). We have A^, = Xn/L{an), 
and with L{l/q) = ^{q)/'p{q), the hypothesis (2.28) ensures L is slowly 
varying and L L. Hence, for any q G (0, oo) we have 

^(-n^) = K^iiM = ^(")(«)%^^ - $(«) 

Han) L{an) 

as n — > oo. By Theorem 3.1, it follows G^"^ converges to H, and by the 
extended Levy-Khintchine representation theorem, this yields F = F va. 
(10.2). This contradicts (10.1) and finishes the proof. □ 



10.2 Self similcir solutions and domains of attraction 

The self-similar solutions are the simplest examples of eternal solutions. All 
self-similar solutions are generated by the s-measures H{x) = Cx^~P with 
p G (0, 1], C > 0, with corresponding Laplace exponents 

m=Cq''^^—^, nQ) = Cq^+P^^^^^. (10.5) 
p pl + P 
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Thus, there is a one-parameter family up to trivial scalings. By (4.8) and 
(5.19), for appropriate C, the Laplace exponent (p = ip{t,q) of the solution 
satisfies 

(f = = — ^ (K = 2), (10.6) 

q = ^p + ^^{e^^) = ^+{e^^Y+P {K = x + y). (10.7) 

The self-similar solutions were described in [16], and can all be captured by 
expressing the associated probability distribution in the form 

F(t,x) = Fp^^{x / dp^^ it)), (10.8) 

for 7 = 0, 1, 2, where the scale factors are 

ap,o{t)=t'/p, a^,i(t) = e*/^ ap,2{t) = \t\~^^^ , (10.9) 

with /3 = p/(l + p), and the probability distributions Fp^j are explicitly 

Fp,(x) = F,,,(x) = if tHp^!r(i+*-w)!^.(io.ii) 

k=l ' ^ 

We now restate the characterization of the domains of attraction of these 
self-similar solutions obtained in [16]. We say a probability measure on E 
is nontrivial if it is not concentrated at the origin. 

Theorem 10.1. Let Ft denote the probability measure associated to a solu- 
tion to Smoluchowski's coagulation equations with K = 2, x + y, or xy. 

(a) Assume there is a rescaling b{t) oo and a nontrivial probability 
measure F on E such that Ft{b{t)x) F{x) at all points of continuity. 
Then there is p € (0,1], and a function L slowly varying at infinity 
such that 

Gto{x)= [ y^+^ut,{dy)r^x^-PL{x), x ^ oo. (10.12) 

J{0,x\ 

(b) Conversely, assume (10.12) holds. Then there is a rescaling b{t) — > oo 
such that 

lim dist {Ft{b{t)dx),Fpj{dx)) = 0. (10.13) 

t— >oo ' 
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Theorem 10.1 illustrates the rigidity of scaling limits. If we insist on the 
existence of a proper limit as t ^ cx) (as opposed to subsequential limits), the 
only possibility is that F{x) = Fp^y{ax) for some p G (0, 1] and a € (0, oo). 
(For degenerate limits, see the remark below.) Theorems 2.5 and 2.9 shed 
more light on this result as they clarify the main hypothesis ((10.12) above) 
and allow us to avoid the use of Karamata's Tauberian theorem in the proof. 

Proof. Let us first prove (a). Suppose there is a (possibly discontinuous) 
rescaling b{t) oo such that limt^oo Ft{b{t)x) = F{x) at all points of 
continuity of F. Then F E ^p, so it is associated to a divergent s- measure 
H. Theorem 2.5 (ii) now implies the convergence of the s-measures 
where 

(x) = ~Xa-^Gto {ax), (10.14) 



ait) 



'm {K = 2), 
6(i)e-* {K = x + y), A(t) 
b{t)\t\ {K = xy), 



t [K = 2). 
e* {K = x + y), (10.15) 
{K = xy). 



As we have seen in the proof of Theorem 2.9, a{t) diverges as t — > Tmax in 
each case. Then by (10.14), the Laplace exponent ifQ for Gt,, satisfies 

XMl/a) ^ Hq) (10.16) 

as i — >■ Tmax) where $ is the Laplace exponent of H. Taking t — > Tmax 
along a sequence tn such that A(i„_)_i)/A(t„) ^ 1, by a fundamental rigidity 
lemma [9, VIII.8.3], we infer that the only possible limits are power-laws, 
meaning ^{q) = cq^ for some /> > 0. Since if is a nontrivial s-measure, we 
must have < /O < 1 and c > 0. Moreover we infer ipQ is regularly varying 
at 0, meaning ipo{q) = qfL{q), where L{aq) / L{q) — > 1 as g — > for every 
a > 0. Note that by (10.16), 

A - caP/Lil/a), c„ = A(t„)99o(l/«(*n)) ^ c. (10.17) 

With tn as described and ctn — a(tn 

), we claim a^+i/a^ ^ 1 as n ^ oo. 
Let a > 1 and suppose a„_|_i/a„ > a for infinitely many n. Then since ipo is 
strictly increasing, along this subsequence we have 

Cn+i A(t„) v3o(l/an+i) , v?o(a"Van) -p ^ 1 
— < -— — ^ ^ a ^ < 1. 



A(i„+l) Cn (foil/ an) (foil/ an) 

But the left-hand side converges to 1. Hence limsupan+i/o^n ^ 1- Similarly 
we deduce liminf a„+i/a„ > 1, establishing the claim. 
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We may now apply the rigidity lemma [9, VIII.8.3] to (10.14) to infer 
that Gtg is regularly varying at oo, meaning (10.12) holds. (The value of p 
must be the same here, due to (10.17) and (10.5).) This proves part (a). 

To prove the converse, we assume that (10.12) holds. Since (10.12) holds 
we may choose increasing rescalings a{t) oo and X{t) such that the s- 
measures G^^\x) = Xa~^Gtf^{ax) converge to H = x^~^. Let b{t) be defined 
by (10.15) for the various kernels. It then follows that Ft{b{t)x) Fp^^{x) 
for every a; > 0. Since the metric is equivalent to weak convergence we also 
have (10.13). □ 

Remark 10.1. A remaining nontrivial possibility discussed in [16] is that 
of a defective limit on E, which we may now take to mean that Ft{b{t)x) — > 
F{x) where F is a probability measure on £^ = [0, oo], with < F(oo~) < 1, 
meaning that gel appears in the limit. If this is the case, then F is an 
element of the full scaling attr actor A, and by Theorem 9.6, the rescaled 
s-measures ^ H, the s-measure associated to F. Moreover, y ^H{dy) 
must have nonzero charge /i^o at oo, and hence $(0+) > 0. This means that 
the in the proof above, the rigidity lemma must yield p = 0, i.e., we must 
have $(g) = c > 0, corresponding to an eternal extended solution consisting 
of a pure dust/gel mixture. By the discussion in [16] (see Remarks 5.4 and 
7.4) a necessary and sufficient condition for this to occur is that 

! y-^GtMv) - L{x), x^oo, (10.18) 

J[x,oo) 

where L is slowly varying at oo. 
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